Prediction of resolvent mode shapes in supersonic turbulent boundary layers

This work applies resolvent analysis to compressible zero-pressure-gradient turbulent boundary layers with freestream Mach numbers between 2 and 4, focusing exclusively on large scale motions in the outer region of the boundary layer. We investigate the effects of Mach number on predicted flow structures, and in particular, look at how such effects may be attributed to changes in mean properties. By leveraging the similarity between the compressible and incompressible resolvent operators, we show that the shape of the streamwise velocity and temperature components of resolvent response modes in the compressible regime can be approximated by applying ideas from wavepacket pseudospectral theory to a simple scalar operator. This gives a means of predicting the shape of resolvent mode components for compressible flows without requiring the singular value decompositions of discretized operators. At a Mach number of 2, we find that accurate results are obtained from this approximation when using the compressible mean velocity profile. At Mach numbers of 3 and 4, the quantitative accuracy of these predictions is improved by also considering a local effective Reynolds number based on the local mean density and viscosity.     

Stokes equations under nonlinear slip boundary conditions coupled with the heat equation: A priori error analysis

In this work, we consider the heat equation coupled with Stokes equations under threshold type boundary condition. The conditions for existence and uniqueness of the weak solution are made clear. Next we formulate the finite element problem, recall the conditions of its solvability, and study its convergence by making use of Babuska–Brezzi's conditions for mixed problems. Third we formulate an Uzawa's type iterative algorithm that separates the fluid from heat conduction, study its feasibility, and convergence. Finally the theoretical findings are validated by numerical simulations.     

竖直恒温面自然对流边界层不稳定性研究

在无扰动、随机式扰动以及正弦式扰动下,通过对竖直恒温面处状态Ra为1.328×10^9、Pr为6.24的自然对流进行模拟,探索了热边界层的不稳定性和共振强化自然对流换热。结果表明:(1)竖直自然对流边界层上游位置的随机式扰动对热边界层的影响主要体现在稳定阶段;(2)该状态下的竖直自然对流边界层的特征频率为15 067,且相比于无扰动状态,频率为15 067的正弦式扰动能在竖直恒温面处提高5.15%的换热量;(3)在竖直自然对流边界层上游位置加入特征频率的正弦式扰动,竖直恒温面处的局部努塞尔数Nu均出现明显波动,且波动随着边界层高度的增加而增大。     

Meshless acoustic analysis using a weakly singular Burton-Miller boundary integral formulation

The Burton-Miller boundary integral formulation is solved by a complex variable boundary element-free method (CVBEFM) for the boundary-only meshless analysis of acoustic problems with arbitrary wavenumbers. To regularize both strongly singular and hypersingular integrals and to avoid the computation of the solid angle and its normal derivative, a weakly singular Burton-Miller formulation is derived by considering the normal derivative of the solid angle and adopting the singularity subtraction procedures. To facilitate the implementation of the CVBEFM and the approximation of gradients of the boundary variables, a stabilized complex variable moving least-square approximation is selected in the meshless discretization procedure. The results show the accuracy and efficiency of the present CVBEFM and reveal that the method can produce satisfactory results for all wavenumbers, even for extremely large wavenumbers such as k = 10 000.     

Numerical investigation of particle lateral migration in straight channel flows using a direct-forcing immersed boundary method

Inertia-induced cross-stream migration has been recently exploited for precise position of particles in confined channel flows. In this work, a three-dimensional finite volume based immersed boundary method has been developed to study the lateral migration and hydrodynamic self-assembly of neutrally-buoyant particles in pressure-driven flows. Simulation results show that, in 2D channel flows, the equilibrium position for a circular particle is closer to the centreline for larger particle Reynolds number due to the increasing flow rate, while in 3D square duct flow, the equilibrium position for a spherical particle is near a face centre and is closer to the wall for larger particle Reynolds number. Self-assembly of a pair of particles is observed in 3D square duct flows but not in 2D channel flows. Mechanisms for the self-assembly are discussed.     

Acoustic radiation force on an elastic cylinder in a Gaussian beam near an impedance boundary

Acoustic radiation force (ARF) is studied by considering an infinite elastic cylinder near an impedance boundary when the cylinder is illuminated by a Gaussian beam. The surrounding fluid is an ideal fluid. Using the method of images and the translation-addition theorem for the cylindrical Bessel function, the resulting sound field including the incident wave, its reflection from the boundary, the scattered wave from the elastic cylinder, and its image are expressed in terms of the cylindrical wave function. Then, we deduce the exact equations of the axial and transverse ARFs. The solutions depend on the cylinder position, cylinder material, beam waist, reflection coefficient, distance from the impedance boundary, and absorption in the cylinder. To analyze the effects of the various factors intuitively, we simulate the radiation force for non-absorbing elastic cylinders made of stainless steel, gold, and beryllium as well as for an absorbing elastic cylinder made of polyethylene, which is a well-known biomedical polymer. The results show that the impedance boundary, cylinder material, absorption in the cylinder, and cylinder position in the Gaussian beam significantly affect the magnitude and direction of the force. Both stable and unstable equilibrium regions are found. Moreover, a larger beam waist broadens the beam domain, corresponding to non-zero axial and transverse ARFs. More importantly, negative ARFs are produced depending on the choice of the various factors. These results are particularly important for designing acoustic manipulation devices operating with Gaussian beams.     

螺栓松动边界下纤维增强复合薄板的振动测试方法

提出了基于预埋压力传感器的量化测试方法,研究了螺栓松动边界对纤维增强复合薄板振动特性的影响。首先,自主设计并开发了带有预埋压力传感器的螺栓松动边界下复合薄板的振动测试系统,并详细介绍了系统各个部件的组成和功能;然后,归纳出一套合理、规范的松动边界下复合薄板的振动测试流程,并对HF10碳纤维/树脂复合薄板进行了实际测试。结果表明:随着螺栓松动程度的不断增加,复合薄板的固有频率逐渐降低,模态振型的节线位置也发生了不同程度的变化,但其阻尼结果呈现先增大后减小的趋势;而共振和非共振响应呈现先减小后增大的趋势。     

Space-time finite element methods stabilized using bubble function spaces

ABSTRACT

In this paper, a stabilized space-time finite element method for solving linear parabolic evolution problems is analyzed. The proposed method is developed on a base of a space-time variational setting, that helps on the simultaneous and unified discretization in space and in time by finite element techniques. Stabilization terms are constructed by means of classical bubble spaces. Stability of the discrete problem with respect to an associated mesh dependent norm is proved, and a priori discretization error estimates are presented. Numerical examples confirm the theoretical estimates.     

Scattering of Tollmien-Schlichting waves by localized roughness in transonic boundary layers

The laminar-turbulent transition in boundary-layer flows is often affected by wall imperfections, because the latter may interact with either the freestream perturbations or the oncoming boundary-layer instability modes, leading to a modification of the accumulation of the normal modes. The present paper particularly focuses on the latter mechanism in a transonic boundary layer, namely, the effect of a two-dimensional(2 D) roughness element on the oncoming Tollmien-Schlichting(T-S) modes when they propagate through the region of the rapid mean-flow distortion induced by the roughness. The wave scattering is analyzed by adapting the local scattering theory developed for subsonic boundary layers(WU, X. S. and DONG, M. A local scattering theory for the effects of isolated roughness on boundary-layer instability and transition: transmission coefficient as an eigenvalue. Journal of Fluid Mechanics, 794, 68–108(2006)) to the transonic regime, and a transmission coefficient is introduced to characterize the effect of the roughness. In the sub-transonic regime, in which the Mach number is close to, but less than, 1, the scattering system reduces to an eigenvalue problem with the transmission coefficient being the eigenvalue; while in the super-transonic regime, in which the Mach number is slightly greater than 1, the scattering system becomes a high-dimensional group of linear equations with the transmission coefficient being solved afterward. In the largeReynolds-number asymptotic theory, the K′arm′an-Guderley parameter is introduced to quantify the effect of the Mach number. A systematical parametric study is carried out,and the dependence of the transmission coefficient on the roughness shape, the frequency of the oncoming mode, and the K′arm′an-Guderley parameter is provided.